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TECHNICAL NOTE 2705 


theory of sufersonic potential flow in turbomachi nes 

By Robert H. Wasserma^ 

SUMMARY 

A generaL method for solving supersonic potential flow problems for 
stationary or rotating coordinate systems is presented. The principal 
attributes of the method are: It can handle flows which cannot be 

treated as two-dimens ional, and a sound theoretical basis gives assur- 
ance of its validity for a class of boundary- value problems. An appli- 
cation to the design of a compressor rotor is made. 


INTRODUCTION 

The fluid flow through a turbomachine is intrinsically three- 
dimensional. This fact must be considered for a full understanding of 
such flow , and in particular for adequate treatment of such problems as 
off-design performance surging and secondary flows due, for example, to 
boundary layers. However, this three-dimensionality is not easily 
accounted for theoretically, even in the idealized case of no viscosity 
or heat transfer. In existing approaches, some sort of two-dimensional 
flow is first considered, such as flow through a cascade (Tyler, refer- 
ence 1), axially symmetrical flow (Marble, reference 2; Goldstein, 
reference 3) , flow in surfaces of revolution (Wu and Brown, refer- 
ences 4 and 5) , or flow over other special surfaces (Stanitz, refer- 
ence 6) . In some cases, this flow serves as a first approximation and 
is modified to give normal variations by Taylor series expansions 
(Reissner, reference 7), by use of Ackeret's two-dimensional vortex-and- 
source method (Meyer, reference 8) or by successive applications of two 
or more different types of two-dimensional flow (Wu, reference 9) . 

The solution of the two- dimensional problems and the extensions to 
three-dimensional flow both generally involve one or more of the following 
numerical techniques: use of formal series; repeated substitution of 

"approximate solutions" into the differential equations; and replacement 
of differential equations by difference equations and subsequent appli- 
cation of relaxation or matrix methods. There is no assurance that 
numerical results obtained by such means correspond to a solution of the 
original three-dimensional problem. 
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The present treatment is a direct attack on the supersonic three- 
dimensional problem; it dispenses with special two-dimensional flcwB. 
Moreover, the mathematical basis of the present treatment due to 
E. F. Titt (reference 10) contains a guarantee that the numerical pro- 
cedure involved shall converge “loc al 1y " to the correct ans wer - 

The fluid flow through a single component of a turbomachine is 
treated herein - either a rotating or a stationary component. Moreover, 
consideration is limited to regions in a component in which the fluid 
flow may be considered inviscid, isentropic, and irrotational. Such 
regions, however, are permitted to be bounded by surfaces on which these 
assumptions are not valid - such as shock surfaces . In such a region, - 
the flow is described mathematically by 
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when the flow space is provided with cylindrical coordinates z,cp,r 
which rotate with the angular velocity cd of the wheel. Note that the 

sound speed a is a function of I, a>, r, ^ and A 

solution of this equation, and thus of a flow problem, is obtained in 
the region when suitable physical conditions are given on the boundary 
of the region. A description of a method of solving equation (l) with 
suitable bo undar y or initial condition will be presented first from a 
general point of view. Then an application of the method to a specific 
flow problem will be made. The work described herein was done at the 
NACA Lewis laboratory. 
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GENERAL THEORY 


For the present purposes, equation (l) is characterized hy the 
statement that it is a single quasi- linear partial dif ferential equation 
of the second order for a single function of three independent variables; 
that is, it is an equation of the form 


11 d 2 $ -12 

a r — — + 2a 


dz2 


a 2 $ 

Sz Sep 


+ 2a. 


13 


a 2 $ 

dz dr 


+ a 


22 


a 2 $ 

Scp2 


+ 2a 


23 


Sep dr 


+ a 


33 d 2 $ 


dr 2 


+ h = 0 


where a^ and h are functions of z, qp, 

(All symbols are defined in the appendix.) 
for solving equation (2) is a so-called characteristic method. The 
central idea of such a method is to solve a differential equation hy 
replacing it hy an equivalent system of differential equations, each 
equation containing derivatives with respect to fewer independent 
variables . This is done hy utilizing characteristic manifolds. The 
problem is thus reduced to solving the equivalent system called the 
system of characteristic equations (Courant and Hilbert, reference 11) . 
In reference 10, a system of characteristic equations for equation (2) 
is obtained by means of characteristic surfaces, and a constructive 
existence theorem for the characteristic equations is presented. 
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Characteristic Surfaces 


On any surface 


z = z(s,w) 
cp = cp(s,w) 
r = r(s,w) 

the following six second-order strip conditions are required: 

a b$ _bh dz a 2 $ acp a 2 $ dr . 
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2 2 2 
d d$ d $> dz d $ dcp d $ dr 

33 3^> dcp dz dw 1 ^qj2 dw 1 dcp dr 33 

(5) 

d d$ d^$ dz d%> dcp d^$ dr 

33 3r - dr dz dw dr dcp dw ^ r 2 33 

(6) 

2 2 2 
d dJ d $ dz , d $ dcp t d $ dr 

33 33 - ^ z 2 ds 1 dz dcp ds dz dr 33 

(7) 

2 2 2 
d d$ d $ dz d $ dcp d $ dr 

33 3<p - dcp dz 33 1 ^cp2 ds 1 dcp dr 33 

(8) 

d d$ d^$ dz d^$ dcp d^$ dr 

33 3r ~ dr dz ds dr dcp ds ^ r 2 33 

(9) 


Since the determinant of the coefficients of the partial derivatives of 
$ of the right-hand side is zero, this is a singular system. If the 
determinant of the coefficients of equation (2) and the coefficients of 
the five linearly independent expressions on the right-hand sides of 
equations (4) through (8) is zero, that is: 
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then the surface given by equation (3) is called a characteristic 
surface . Clearly, whether a surface is characteristic depends on the 
solution $(z,qp,r). Ihe left-hand side of equation (10) can be expanded 
by Laplace's expansion by minors of the last two rows to give 


Al2 (If) + A 14 If If - (A 15 + A 23>|f|f- 
A 24 (g) 2 4 u 25 ♦*4>g£-* 56 (I ) 2 = 0 M 

where the Aj_j are minors obtained from the first four rows, with the 
i^k and col umns omitted. 

Just as the two- dimens ional method requires the existence in the 
plane of characteristic curves (see- ref erence 12), this method requires 
the existence in space of characteristic surfaces. Thus, the class of 
solutions of equation (3-1) is examined. In this class, there are those 
characteristic surfaces which can be parametrized by two space coordi- 
nates. For concreteness and because they are suitable for later appli- 
cation, the following two cases are considered: (l) s = cp and w = r; 

and (2) s = z and w = r. In case (l) , equation (11) becomes 


A 32 (§f) + A 14 || - ^24 = 0 ^ a ) 

Cp ( C^Z C^P C^I*\ ( \ 

where A-^j are obtained from A^j by putting ygfi* ~ l 0^ 11 

and taking the elements of the first row of equation (10) as functions 
of cp and r. In case (2), equation (3.1) becomes 

.z .z Sep .z [ Sep\ 2 

A 12 + A 14 37 - A 24 ) = 0 ( 1213 ) 

z / Sz Sep Sr\ f A Sep 

where Aij are obtained from A-y by putting \^, ^ j = 1<J, 1 j 

and taking , the elements of the first row of equation (10) as functions of 
z and r . 


Since equation (12a) is a first-order equation in two independent 
variables, it may be solved with a given noncharacteristic initial curve 
z = f^(t) , cp = fg(t), r = ^3(^) having first derivatives such that 
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df-3 


Bz 


Bz 


^ ^ 0, such that initial values of ^ and ^ satisfy equation (12a) 

on the curve, «nd such that 


“3 

dt 


Bz ^2 


Bz ^3 


dt + "Sr dt 


(13) 


(see reference 11, p. 63 ff.), and, similarly for equation (l2h) . Then 

Bz Bcp 

for these two cases, respectively, the initial values of ^ and ^ 
must satisfy 


/ \2 

.* ( Bz\ . * Bz .* _ 

^12 \m>) + A 14 " -*24 = 0 


(14a) 


♦ * Sep N * / Sep\ _ 

a 12 + A 14 -*24 l^j = 0 


(14b) 


where A-y are obtained from A^j hy putting 

^ §) takiD s the elements of the fir8t row of 

equation (10) as functions of t. Equation (14a) has two solutions for 
Bz 




if 


A 14? + 4 A 12 ^24 > 0 


(15) 


4 2/0 

Equation (l4h) has two solutions for ^ if 


(16a) 


a 14 2 + 4 a 12 ^24 > 0 


A 24 / 0 


(15) 

(16b) 


Thus, u nd er either of these pairs of conditions there are two character- 
istic surfaces through the initial curve. 
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Recall that in the two-dimensional case, there is a criterion for 
the existence of two one-parameter families of characteristic curves in 
the plane. These are taken to form a new curvilinear coor dinat e net - 
the characteristic coordinate net. The characteristic equations are 
then obtained as equations in the characteristic coor dinat es or param- 
eters '(reference 12) . In three dimensions, the criterion is for the 
existence of two characteristic surfaces through each curve of a set of 
curves in space. If appropriate families are selected from this collec- 
tion of surfaces, the procedure is somewhat analogous, l&king such a 
selection requires specifying an initial value problem. 


Initial Value Problem 

Since, for example, $ = constant satisfies the partial differential 
equation (1) in an arbitrary region, it is clear that a solution of 
equation (2) satisfying further conditions is wanted. Thus, in addition 
to the condition that a function $(r,cp, z) satisfies equation (2) in a 
given region, further conditions are required of the farm: on a speci- 

fied part 2 of the surface bounding the given region, $ as well as 

its first derivatives ^ are to take on a priori given values. 

These may be given arbitrarily up to certain limiting conditions. 

Let an initial surface 2 be given. It is assumed to be not char- 
acteristic ; that is, it does not satisfy equation (ll) . A one-parameter 
family of curves is selected which simply cover 2 and which are nowhere 
tangent to r = constant curves; each curve may therefore be parametrized 
by r. Under the conditions (15) and (16a) or (15) and (16b), there is 
a pair of characteristic surfaces passing through each curve, and these 
surfaces plus r = constant surfaces may be taken as coordinate surfaces 
of a new coordinate net, since the Jacobian of the transformation 


r (17) 


J 


z = z(u,v,r) 
cp = cp(u,v,r) 
r = r 


is not zero on 2 (reference 10) . 2 may be given by 
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z = z(s,r) N 

q> = cp( B ,r) y (18) 

r = r 

where s Is now the parameter of the family selected, and the character- 
istic surfaces are 


z = z(u,v Q ,r) ^ 
cp = cp(u,v 0 ,r) l 


and 


r = r 


J 


z = z(u Q ,v,r) 
V = ^o^r) 


r = r 

For convenience, put u = s and v 

* \Tz 


J 

-u on 2. 


(19) 


( 20 ) 


The limiting conditions on 2 other than those already mentioned 
are; the functions defining 2 and their first derivatives with 
respect to s are partially analytic functions' with respect to w. 

This is to he true of all the functions defined on 2 (cf. reference 10). 
The remaining limiting conditions on the initial functions are the strip 

S$ S$ 

conditions. On any surface, in particular on Z, must 

satisfy the two first-order strip conditions 

S - Sz 3$ Sep S$> Sr 
+ + 3?3i‘ 


( 21 ) 
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S , _ S$ Sz S$ Sep S$ Sr 
Sw - Sz Sw + S^p Sw + Sr Sw 


( 22 ) 


Characteristic Equations 

In terms of the parametrizations of equations (19) and (20), the 
characteristic surfaces are obtained from equation (11) in the forms 


u ( Sz") 2 u Sz Sep u /Se^ 2 ^ . 

*■12 + A 14= Bu Su " A 24 [^) = 0 ( 23 ) 

■v ( Sz\ 2 v Sz Sep v / Sep \ 2 

A 12 \^F/ + A 14 SFF ' A 24 \&) = 0 (24=) 

u v Sz Sep 

Note that and A^j are the same functions of a , the 

distinct notation indicating that here in the first case a^ are con- 
sidered functions of u,r and in the second case a^ are considered 
functions of v,r so that in general the common notation will 

suffice. 


a 12 / 0 implies that / 0 on 2. Jet the root determined by 

. . +/ fi 14 + 4Q L2 ^24 

equation (23) considered as a quadratic he = 2^2 

and let the root determined hy equation (24) he 
_Q L4 ~/ Q 14 + 4Q L2 Q 24 


p 2 = 


14 

2^L2 


Then equations (23) and (24) become: 
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Similarly, if Ag4 j- 0 then equations 


(23) and (24) became 


S*P _ Sz 
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(25b) 
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Sep Sz 

SE “ T 2 ^ 


0 


(26b) 


where t,- = — . 

1 P± 


On any surface the 12 third- order strip conditions must also hold. 
They are arranged as follows, .where e takes on the values z,<P,r 
successively and subscript notation is used for derivatives with respect 
to z,cp,r: 
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S . , Sz . Sep . Sr 
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S . . Sz * Sep , Sr 

SE V s ®zepe SE + 9 W*SS + ^repe SE 

S . . Sz . Sep . Sr 

SE $ r£= $zre SE + $ cpr£ SS + **ve SE 


In addition, the third derivatives must satisfy 


(27) 


alI ^ Z Z£ + ZaZ1 %z € + 2a31 $rze + &2 %<pe + 2a ^ 32 $ r <p € + a 33 $ rre + a^ zz + 

2a 21 $ C p Z + 2a 31 $ rz + a? 2 ^+ 2a 3 £ 2 $ rcp + a 33 ^ + b £ = 0 (28) 

which are obtained by differentiating equation (2) . Let 
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B(e) = A 11 ^ + 2a 21 $ + 2a 51 $ + a 22 $ + 2a 32 <f> + a 33 $ + b £ 

£ zz £ cpz £ rz £ qxp £ rep £ rr 

'■ ' (2d) 

so that in equations (27) and (28) there are for each £ seven linear 
nonhomogene ous equations in six unkn owns which may he written as: 
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Wow if the surface is a characteristic surface, the matrix of equa- 
tion (30) is of rank 5 or less ( c f. equation (10)) and thus, as consist- 
ency requires, so is the augmented matrix. This in turn requires that 
the two sixth-order -min ors formed from the augmented matrix by deleting 
the last row and second column and the first row and second column be 
zero. In particular, on the characteristic surfaces given by equa- 
tions (19) and (20), these minors upon expansion become: 
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on t = constant, 
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on u = constant. 

By means of equations (25a) and. (26a), equations (3l) and (32) 
expressed in terms of P-^ and p^: 


Q 24 5 “® z €-%2 Pi |l V + E( £ ) || = 0 
Q 24 |f ®ze - Rl 2 p 2 ^7 V + E(e) §|=0 


(32) 

(33) 

(34) 

may "be 

(35a) 

(36a) 


By means of equations (25b) and (26b), equations (31) and (32.) may Be 
expressed in terms of T and Tg: 
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Now consider equations (25b) and (26b); equations (2l), and (7) 
to (9) with s = u; equations (35b) and (36b) with z for £ ; equa- 

tion (35b) with (p and r for £ ; and equation (33) with cp and r 
for £. These are 12 partial differential equations for 12 functions 

(y«x) = (z,^ ^^ z j^cp^ r ^ Z z^ziCp-» l ^zr- , ^cpcpJ^cptr-’^ > rr) 

in (u,v, r) space. They are first-order nonlinear equations, which, 
however, are linear and homogenous in derivatives with respect to u 
and v. They have the form 
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> (37) 
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where and are functions of y a and 


In the (u,v,r) space the initial surf ace - 2 "becomes u + v = 0. 
The values of the functions y a on 2 are (y a )(°) = y a (u, -u,r) = 
y a (^2 u,r), where y a are initial values on 2 as fractions of s,r. 

_ . _ ^o. „ (*«)«» *y a 

The values of ^ — on z. are - 


c3r 


5F7 


5F** 


^y a 


^yaV 0) 

The values and ' oi 

obtained by solving equation (37) on the initial surface with 


of «s — and ^ — on 
dv 
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a i r k\- 0) 

[5Tj = ^ ST 


k = 1, 2, 3, . . . , 12 
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(38) 


Equations (38) come from y a (s,r) = y a (u(s), v(s), r) and the relations 

d^y a d^y a 

s=/2u=-/2v on the initial surface. Finally, and ^ 


on 2 are, respectively 


p 2 y a Z (0) _ bfiy a ) 

3 \dr du / 5r \du / 


(0) (* 2 y*\ i0) 

“ a \srsv - 
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The 12 equations (37) with these initial functions on u + v = 0 
constitute an ini tial value problem. It follows from the proof of refer- 
ence 10 for the general case that this problem is equivalent to the 
initial value problem for equation (2) . A solution of this problem thus 
leads imm ediately to a solution of the initial value problem for equa- 
tion (2) . Uniqueness is also preserved. 


Solution of Equivalent Problem 


Equations (26b) and (36b) with e = z are differentiated with 
respect to u, and the others are differentiated with respect to v so 
that there are obt ain ed 12 equations linear in 12 second- order uv 
derivatives and nonhomogene ous , 
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in which the elements 

of the empty rectangles are 

all 

zero. 

Briefly, 
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The augmented column Gj_ of equation (39) is 
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Solution of this system yields a set of equations of the form 


d 2 y 
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du dv a 


a = 1, , 


12 


(40) 
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where 


/ *3 Sy 3. S2y 3 * y 3 

fa - fa 1 yp, ^T> 57 -, 5TSP ^T5F 


a 2 


w 

since Tg = 0 when =* then / ^Dsa ^ u "^ v = Gg is satisfied 


d 2 ; 


a 


ye 


identically and may he chosen arbitrarily and the 

equations solved for the remaining second derivatives. To solve the 
transformed initial value problem, reference 10 uses an extension of the 
Picard method. 

If equation (40) is integrated with respect to v from the initial 
surface u + v = 0 to an arbitrary point u,v,r (see fig. 1), 


dy. 




ST (u,v,r) - (u,-u,r) = f a dv 


(41) 


Then integration again with respect to u from u = -v to any value of 
u gives y a as 

r Sy a 

y a Kv,r) = y a (-v,v,r) + / (u,v,r) du 

\J -V o 


= ya(-^ T ^ r ) + 




T K-u,r) + r ^ 

J -v L vJ-U 


dv 


du 


(42) 


The quantities y a (-v,v,r) and — — (u,-u,r) are known, but f a depends 

chi 

on values of y a at points- not on u + v = 0; a zeroth approximation 
r (0) 


u a 


is therefore chosen for f_ and a method of successive approxi- 


mations is employed. Let f 


( 0 ) 


a 


be f a with the values of its argu- 


ments assuming their initial values at corresponding v. This gives new 

values for (u,v,r) and in turn new values (y a ) ^ 

y a (u,v,r) . The relation 


for 
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dy. 




PU 


"a t N “a , N 

Sv" (u,v,r) = (-v,v,r) + 


f a du 


(43) 


-v 


^aV 1 ) , (i) /^NC 1 ) /^y \(!) 

similarly glares . The functions (yj ^\ i '5uT/ ' \¥7 

and their derivatives "with respect to r are substituted back into 
equations (42) and (43) and, the process is repeated. In general. 


dy a \0>+i) 
55“. 


SI 


pv 


fa (W> IV 


(44) 


-u 


ZyMv+1) /dy\(0) 

377 = V377 + 


pu 


du ' 


(45) 


and 


IU p y „V 0+1 > 


. (v+1) (0) P /^ y ct) 

y <* - y a w + ['srJ 

VJ -V 


du 1 


(46) 


where v = 1, 2, . . . 

The y a ^ converge to the solution of the transformed prohlem. 


In the application which is now to be considered, the design of a 
rotor, AjJ^_ / 0 does not alwayB hold, so that the equivalent problem must 

also be expressed in terms of equations (25a) , (26a) ; equations (21) , and 
(7) to (9) with s = uj equation (35a) with z and cp for e; equa- 
tion (36a) with cp for e; equation (35a) with r for ej and equa- 
tion (33) with z and, r for e. The use in this application of two 
different formulations of the equivalent problem can be avoided by 
putting w = cp or w = z instead of w = r when characteristic sur- 
faces are selected. Then, the physical condition which must be satisfied 

is w 2 + wjp> a 2 , which may be assumed to hold everywhere. 
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ROTOR DESIGN 

Shock Surface E^ and Upstream Suction Surface Eg 

As an application of the general theory, the flow of air through a 
rotor of an axial compressor is considered. Attention is focused on a 
single rotor passage and on the air flowing into, through, and out of 
this passage. As the air upstream of the "blade passage flows into the 
"blade passage, the "blades affect it in one of the following general ways. 
In the idealized case of thin "blades, the upstream velocity can "be arbi- 
trarily prescribed subject only to the condition that the components be 
uniform in the tangential direction, and then if the design point is 
chosen so that a blade is tangent to the relative stream surface at its 

, 1 £cp cta> 

leading edge (that is, if the blade satisfies — w ^ = w r + w z 

along its leading edge) , then the air will flow smoothly into the 
passage. On the other hand, if the blades are not thin (nor have cusp- 
shaped leading edges) then for no steady upstream velocity that may be 
specified can both the pressure and suction surfaces satisfy the condi- 
tion, so it would seem that there would have to be a discontinuity in 
velocity as the air enters the blade passage. An indication of what 
actually happens for a blade with arbitrary leading edge is given by the 
following result (cf . Jones, reference 13) for approximately straight 
leading edges with approximately constant wedge angles. 

(a) If the component of the upstream relative velocity in the plane 
normal to the leading edge is supersonic, then in all other respects the 
upstream velocity may be arbitrarily prescribed and a discontinuity, or 
shock, generally will occur. 

(b) If the component of the upstream relative velocity in the plane 
normal to the leading edge is subsonic, then this velocity component 
makes zero angle with the blade in this plane and the air passes smoothly 
into the passage. In general, the condition on the velocity . in (a) 
permits specification of upstream velocities within a considerable range 
including steady, uniform upstream velocities . In (b) the rotor effects 
an adjustment of the velocity upstream making an a priori prescription 

of the velocity upstream unrealistic. 

In view of the foregoing considerations, it will be assumed in the 
sequel that the upstream velocity components are prescribed according 
to (a) . If the wedge angle (at each radius) is small, enough and if at 
design point the suction surface at its leading edge is tangent to the 
upstream relative flow, then there will be 'a shock surf ace Z ^ coming 

off the pressure surface attached to its leading edge © and going 
across the passage and downstream toward the suction surface with leading 
edge (2) (fig. 2) . 
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The shock surface is described by 


r 2 

f'vr \ 2 

(— ) -X 

Cpl + r 2 1 

( 


L W J 


_\ a/ _ 


„2 _ 2 w z w r ~ „ 

q>r + 2r 


v /v . 

2r -^V 2 - 2r <P r + (-£) - X =0 


(47) 


with a given leading edge. In equation (47) , subscript notation is used 
for derivatives of cp; X is a measure of the entropy change across 
and the coefficients are evaluated upstream of 2^. In order to 

obtain the differential equation for 2-^ first note that 


cos = 


cp r -w r - - w + °PzV z 




~ v q> + ^Vr 
gw 


w r + w cp + v z 


(cf. fig. 3). Or in terms of p-^. 


(«p w - w + rcp z w ) 2 
?? ” =Sln Pl 


(47a) 


But by means of the conservation principles for shocks and the Erandtl 
relation, sin 2 . P-j_ is in turn expressible in terms of pressure ratio 
according to 




-§ sin 2 p! - p 2 
a l 


(47b) 


(see reference 12.) 

Finally, the Rankine-Hugoniot relation gives 

pJ 


p 2 - K 


p 2 E 




Pi/ 


1 
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where K = e ~^ S / C P. That is, 
change across the shock. 


P2 

— may he expressed in terms of entropy 



(47c) 


when sin d 0-. and. — are eliminated among equations (47a) , (47b) , and 

Pl f / K \ . ..2 

(47c) , equation (47) is obtained with — 




r(i - ii ' ) 


It is assumed that Z 2 intersects the suction surface along a curve 
C 22 - Then the blade passage is split into a region upstream of 2^ 
and a region downstream of Again, for simplicity, let the suction 

surface between the leading edge (2) and the curve c-^ be constructed 

to coincide with a stream surface of the upstream velocity. This part 
S 2 of the suction surface then is described by 


VP Z + VRr " | “Wcp = 0 ( 4 = 0 ) 

with a given leading edge, and the intersection of Z 2 and 2g deter-' 
mines Cj^ • Thus, the effect of the rotor is confined; to the region 
downstream of Zq. 

The entropy is now assumed to be uniform upstream and downstream of 
Iq_, and its change is constant across £ 2 * Then if the flow is irrota- 

tional upstream of it will be irrotational downstream of 2^. As 

a result, equation (l) will hold downstream of The downstream 

irrotationality under the given assumptions may be demonstrated as 
follows: On either side of 22 (reference 9) 


- WX(VXV) = - VI + WS 
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By hypothesis, VI = 0 upstream of 2^ Since I does not change in 
cross ing a shock arid entropy is constant downstream "by hypothesis, 

VI = 0 downstream of (use T ^ reference 9) . Thus downstream 

of 2-^ WX(7XV) =0. If VXV / 0, then a vortex line must follow a 
relative stre amlin e. Since it cannot terminate, it must intersect 2-^ 
and so V* V has a component normal to 2^ at 2^. However, this Implies 

that the vorticity upstream is nonzero aLnce the normal component will 
not change across the shock. Thus, VXV / 0 leads to a contradiction. 

Note that among the consequences of the assumptions X and I are 
constant. For simplicity it will he assumed that H = constant and 
VqT = constant upstream. Moreover, if w r = 0 upstream of 2^, then 

v z = constant and equation (47) is of the form: 

Q 

(C^r 2 + C 2 )Cp 2 + (C 3 r 2 + C 4 )cp2 + (C 5 r 2 + Cg) cp z + C ? r 2 + C Q + ~~ = 0 


' (49) 

Kg 

since a 2 = All C i and are constants. When equa- 

tion (49) is subjected to a Legendre transformation 


F(£,t]) + c P(z,r) = z£ + rrj 


, Bcp 


Bcp 
n = 3r 


z = 


dF 


the following expression is obtained: 


dF 

r = ^ 


(c^ 2 + C 3 T) 2 + c 5 i + c 7 )^ 4 + (C 2 £ 2 + C 411 2 + Cg£ + c 8 )(^J + c 9 = 0 


( 50 ) 
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so that 


F(£,il) =Jit(5,ti) <in + L(5) . (51) 

If the velocity on the upstream side of the shock is 
(v 2 i, w r -j_) and on the downstream side (w z2 , w r2 ) , then 

from the shock relations , 


2 1 - p 2 l-p 2 

w z2 = rCp z a l r<p z w zl - wqj! + tfp r Tr rl + a 2 


-2 


CT l w zl + *Pz v «fiL - ^zFVrl 


•\ 


>2 = ' a l rijy, rl * 72 [ r ^l + °lm + 


2 1 - u. 1 

w r2 = rep^ ^ 


1- 


zl “ w cpi + rcp r w rl CT 


#[- 


^z^z^zl + “Pr% + a 3 w rl 


(52) 

where a-^ is the upstream sonic velocity and a ± are functions of 
rcp r , rcp z related to <7. This is shown as follows. 

If the velocity is decomposed in the direction normal to and 

in the two parametric directions on Z-^, the components beiDg denoted 
by Wjp W TZ , respectively, then on either side of Z^ there is 

the relation (fig. 2) : 


TT „11 .12 .13 ^ 

% = A w r + A Wq, + A w z 

tt .21 .22 .23 . 

W TZ = A w r + A Wqj + A w z > 


tt ,31 ft 32 ,33 

W TR = A . w r + A w cp + w z J 


(52a) 
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■where the matrix of coefficients is 


f 

- a 

rcp z a ' 

a i a i 

«P r a,. 

-rep rep a, 
r ^ z 1 



2 

a l CT 3 


Equation (52a) comes from 


% = 


rCp w - wv + rep v 
z z cp r r 


W TZ - ^ w r “ % ~Er)°l 


Wmo = lW„ - Wl 


where Wjj —— is the projection of ^ in the radial direction, and 
Wpj is the projection of Wg in the axial direction. 


2 2 ,2 

a l = a “ (r%) 


a| = a 2 - (rcp r ) 2 

W Tz/°l is projection of Wrp Z in the radial direction and 

is the projection of Wr^ in the axial direction. Since the deter minant 

of the matrix of coefficients of equation (52a) is not zero, 

n ip "\ 

w r = B W N + B W TZ + B % 


Wq) = B % + B % z + B 


W z = B 31 % + B 32 + B 33 Wrj^ 


(52b) 
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where the matrix of coefficients is 



Inasmuch as equations (52a) and. (52b) hold on both sides of the shock, 
the subscripts 2 may be placed on velocity components in equation (52b) 
and the subscripts 1 on velocity .components in equation (52a) . 

Furthermore, note that W TZ1 = W TZ2 , W N1 W R2 = cl - 

(Erandtl relation, reference 12) and c 2 = p 2 (2h^ + W?) . Substituting 

equations (52a) in (52b) gives equation (52) . 

Downstream Suction Surface 2^ 

As mentioned in the section GEHERAL THEORY, it is necessary to know the 
first derivatives of $, which is essentially the velocity in this case, 
on a surface E in the rotor passage. This surface is taken to be one 
obtained from 2^ and 2^, the part of the suction surface downstream 

of Cj^ • It will be seen that on such a surface the conditions required 
in GEHERAL THE ORY can be satisfied. Thus, conditions on E 3 must now 

be prescribed to completely determine its shape and the velocity on it. 

This may be done in many ways. How it is done best, the factors involved, 
and the limitations imposed depend on performance and constructibility 
requirements of the rotor. If, for the moment, it is assumed that the 
shape of E 3 and the velocity on it have been completely determined, 

then according to the GEHERAL THEORY, the velocity (w r , Wq>, w z ) can be 

determined at every point in the rotor passage. Thus, in particular, 
the pressure surface is obtained as the stream surface through the 
leading edge '!) j the velocity distribution on it, and the velocity dis- 
tribution downstream of the rotor are also obtained. On the pressure 
surface there are certain constructibility requirements, and on the 
velocity distributions, certain performance requirements. An obvious 
requirement is that the pressure surface should intersect the suction 
surface in a curve C34. which may be taken as the trailing edge of the 

rotor blade. Another requirement is that the velocity distribution on 
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the pressure surface (as well as on the suction surface) not give rise 
to shocks and he fairly "smooth" for high efficiency. Further, uniform 
work input or free vortex flow downstream of the rotor is desirable. In 
order to meet these and other requirements, there are the choices in the 
selection of conditions on 2 ^ plus the choice of the as yet undetermined 

constant, X. The following procedure is -given as an example. 

For strength and simplicity of construction, a good specification 
would seem to he that 2^ have straight-line elements slop ing slightly 

toward the radial direction in z = constant planes (fig. 1) . Thus 
2 3 : cp= cp(z,r) would be given illicitly hy 


sto Qa(z) - q>] . - K*) m(s) (ss) 

with the stipulation that m(z) - cp he small anfl positive. Here 2(z), 
the intercept on the line cp = and m( z) , the angle with the line 

C P= 0, are determined in the interval ( z Zg) in terms of cp evalu- 
ated on Cj£ hy 

tan [m(z) - cp] = rq> r (z) . (54) 

which is obtained hy differentiating equation (53) with respect to r 
and 


sin [m(z) - <p(z)J = - X &_£°» (55) 

r(z) 

where c jg is given hy 9(z) and r(z). In general, ■ cp z and cp r are 
determined on cjj> hy 


dcp _ Sep Sq? dr 

dz ~ c3z "Sr dz 


1 Sep 

r w cp = + 3? 


the latter condition being required in order to avoid reflected shock 
along c-j_2 . 
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Beyond Zgj l(z) and m(z) may tie chosen to give the 'suction sur- 
face the desired tur ning . These functions and their first derivatives 
are required to he continuous at z = Zg as well as beyond. Thus 

d , *— P = g(z,r) may be prescribed, where g(z,r). is determined from pres- 
et z 2 

sent exper ime ntal or theoretical knowledge of maximum blade turning 
before stalling. Differentiating equation (53) twice with respect to z 

and eliminating ^ nnri cp yield the relation 

t] r [Z"(z), Z’( z), l(z) , m" ( z) , m’(z), m(z)] = 0 (56) 

Let 

3 2 

m(z) = Az + Bz + Cz + D 

and from the boundary condition on m( z) and m’ ( z) at z = Zg and a 
downstream curve z = z d , the coefficients A,B,C,D are obtained. From 
equation (56), l(z) is then obtained with boundary conditions at z = Zg. 

Tn order to obtain the velocity components on 2 3 , it 1 b first noted 
that on 2 3 they must satisfy 


> = rcp z w z + rcp r w r (57) 

- - 2omcp z = 0 (58) 

Equation (58) i 6 a necessary and sufficient condition that there exists 
a function $>(z,r) on Z 3 which satisfies the first-order strip con- 
ditions. For 2-^ such a condition follows automatically from the 
assumption of upstream irrotationality. 

If the magnitude of the (relative) velocity on is then pres- 

cribed according to what would seem to be a desirable blade loading, then 
eliminating w r and w z among equations (57), (58), and 

W 2 = -w^ + -^ + w^ yields an equation for w^ - of the form 


dw„ dw ro 


rcp„ 
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&l( r > z >V W ) ^ w cp + 8 2( r * z * v qy w ) If w cp + S 3( r ^ z ^qv w ) 


0 

(59) 


In this equation. 


= 2 r ^ ) r <p z ° 2v z 
6 2 ( z^w^W) = 2r 2 cp r cp z o 2 w. 


r~z "r 


6 3 ( z , r ,Vcp,W) = 21^^02^ + 2rcp r rcp z o2w r + (L<P Z v 2 + 


(Uz ^ + rc Rr sr) ^ + (*& sr) - 2rc ?z {[(^) 2 c p z + ^r] 


+ 


(rcp r ) 2 2mrcp z |w a + 2r<^.| (rcp z ) 2 % - + (r<p z ) 3 2o> Jw r 

*Wz ^ (rcpj.) 2 fl 2 - rc^ ^ (rcp z ) 2 W 2 


a| = a 2 - 1 


The -w z and w r appearing in the expressions on the right-hand sides 
are eliminated, respectively, by 

a 2 w z " 2rC Pz v q? r z + CT l4 - (^r) 2 ^ = 0 ( 60 ) 

°z4 - 2z< Pr w cp w r + - (^z) 2 W 2 = 0 (61) 

Alternatively, ^(or vg) may be prescribed on S 3 according to 
a desirable loading. .Then conditions (57) and (58) lead to 
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^z av z ^z fozz ^z^rz ) 
cp r Sz~ 3r~ \cp r <p 2 ) v z 


= -<*> 


drv^ 

dr 


p 

cor -rv. 



rz 


(62) 


Equations 
on C-J 2 - 
z = z d . 


(59) and (62) are solved with given and w z , respectively. 

In particular, a free vortex distribution may be prescribed at 


Validity of Conditions Required in General Theory 

In order to apply the general theory, condition (15) in particular 
must be verified. When the coefficients of equation (2) are taken to be 
those of equation (1) , then 


* 2r' P . . . . p 

^14 = ~ [(r'Wjp - rcp’w r )(r'w z - z’v r ) - a z’rcp’ 


A 12 = “ ]j[( r ' w cp “ rc P , ' w 'r) 2 " &2 (( r ’) 2 + (^’) 2 ) 


4 & = r’[ (r*v z - z'w r ) 2 - a 2 ((r’) 2 + (z*) 2 ) 

When the indicated algebra is carried out, the following expression is 
obtained: 


2 (’-n^ 4 r 


T^W 2 - a 2 ) - (T*W) 2 ] 
where T = (z', rep', r'). Thus, condition (15) becomes 


.*2 4a (r') 

a 14 + 4A 12 A 24 = ~2 

r 


^(W 2 - a 2 ) - (T-W) 2 > 0 
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which reduces to 


cos 2 (T,W) < cos 2 9 (63) 

Since 0 < 9 < jt/2, the angle between T and W must he greater than 
0 , or the curve parametrized by t through a point on ( 2 -^, 2 ^) must 

be outside the Mach cone with vertex at that point. This cond i tion can 
always be satisfied by properly choosing the t curves. If, as in the 
formulation of the rotor problem, - these curves are chosen to be 
z = constant curves, then the application of other conditions such as 
assuming radial leading edges may prevent difficulties. 

If conditions (16a) and (16b) are examined, it is seen that A 12 
and A^ are quadratic forms in rep' ,r' and z',r', respectively. For 

equation (l) , they are definite quadratic forms if w^ + w 2 > a 2 and 

w 2 + w 2 > a 2 , respectively. In the application developed herein, it is 
z r 

as sximed that at least one of these inequalities holds at every point of 
(E^,2g) . Roughly, far enough upstream in the rotor w^ > a and far 

enough downstream in the rotor w z > a. 

Furthermore, it must be shown that the surface (2^,2 3 ) is nowhere 

tangent to a characteristic surface. A sufficient condition for this is 
that the Mach cone with vertex at each point of ( 2 ^, 23 ) cut (2^,2 3 ) , 

because if ( 2 ^, 23 ) were tangent to a characteristic surface at some 

point. It would have to be tangent to the Mach cone at that point. This 
condition is verified by showing that at each point the angle Pg 

between the velocity vector and, its component tangent to ( 2 ^, 23 ) is 
less than the Mach angle. This is obvious for 2 J 3 . For 2-|_ note that 
this condition is equivalent to 


^2 2 

4 Bin 2 P 2 < 1 
®2 


and that the relation 


fl 

P 2 


r(i - n 2 ) 


-| sin 2 p 2 - p 2 

a 2 


(64) 


(64a) 



30 


NACA TN 2705 


analogous to equation "(47b) ■with 


— < 1 becomes equation (64) . 
•^2 


Finally, in order to apply the general theory, the initial surface 
2 has to be parti ally analytic with respect to r in every z interval. 
Clearly, (Z-^Zj) does not satisfy this requirement in intervals con- 
taining points of the curve of intersection of 2^ and 2^. Conse- 
quently, 2 id constructed from (Z-j^Zg) by replacing (2-^,2j) between 
Zj_ and z 2 (cf . fig. 1) -by a surface given by a third degree polynomial 

in z with coefficients functions of r such that it, as well as its 

first derivative with respect to z, matches the original surface on 

z and z„. The same thing is done for the function $. It is 
J- O 

reasonable to assume that the preceding conditions are still valid. 


Integrands for Equations -(44) and (45) 

In solving for the uv derivatives from the linear equations (39) 

^ ^ ^L2 

it is clear that f a is obtained in terms of Q^, Tp, t 2 , -§^~, 

^- 2 E(e), . These must be expressed 

3 3 3 Bv 3 du ' ' ' 3 ov 3 ou 

By 5y a By a By a 

in terms of y a , and derivatives of y a , and with 

respect to r for the integration. 

Evaluating the a 1 ^ as the coefficients in equation (l) gives 


^"7 


/ 1 , \2 
( - cur] 


’) 2 - &2 + 1(1 %- H " ( $2 " &2 ) (i 


(65) 


n 24 = *! - a 2 . - 2 i r *z |f 


+ & 


Bz 

Br 


2 , 


( 66 ) 




Bz, 

Br 


+ $q, - ^ ) $ z - id 


, , Bz - Bcp 

t\t~V - **) ^ BF ■ $ r^z 


( 67 ) 


2400 



NACA TN 2705 


31 


On 2, It vlll be more convenient to have these functions in terms of 

dz 

velocity components. Since now ^ = 0 they become 


- r2 ^L2 = [ w q> " *Pr w r] 2 “ * [j- + far)* 


(65a) 


0 o . = w 2 - a 2 
24 z 


(66a) 


= v z[ v cp “ v r( rC Pr|] 


(67a) 


In order to get the expressions needed for T-^ and Tg, these equations 
are substituted into 


Ql 4 9l4 "V^14 + 4$ \2^24 

T 1 = 2^4 > T 2 = 22 24 


Further, the derivatives of equations (65) and (66) with respect to 

^ T 1 ^ T 2 

u and v are required) for and , these, as well as the 

derivatives of equation (67) , are substituted into 

^24 ^ ^14 - ^14 ^ ^24 

37~ = 


^ 14 ^ ^24 ~ ^24^14^ Ql4 ~ ^24 Q 12 + ^12^4 ^24 

+ 4Q lp24= 


2 d 


^ 24 ^14 “%4 ^24 


2 

= 




+ 


^14 gjj ^24 - ^4%4 ^ ^14 - 2^24 %2 + 2 Ql2^24 Jj^24 


2^4^ Q 2 14: + 40 L2 ^ 4 
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For the E{£), 


E(z) = tf(z) + \_Z%.\ - ($| - a2) |f] $ zz + 

[2 ®cp- a») $r “ (®? " *) M] 3F (® * " &2 ) I? $ zr ( 68 ) 

E(cp) = B(cp) + [~2$ r $ z " (®| " &2 ) ^ $zep + 

[2 ®q,- ®r) 3> r - ($? - a 2 ) $cpcp+ (®? - a 2 ) ^ <£cpr (69) 


E(r) = B(r) + [2$ r $ z - ($ 2 - a 2 ) |f] ^ V + 

[2 V ®) ~ (f? - &2 ) Iffe V + ( $ ? ' &2 ) I? ^ (?0) 

where B(z), B(cp), and B(r) are obtained from equation (29) in which 
a|^ and b e are given by the following: 




“z^zz 


- 2aa z 

a 12 

a z 

*?\u 

II 

r r V r 

- $ $ 
z rz 

a 22 

a z 

II 



f zz 


2 aa z 


z r r r \ r 


= 2t Z 4, Z<(.- 2aa cp 


rz 


a | 3 = 2$ r %. z - 2aa 2 


a cp “ ^r^zcp + ^z^rcp 

23 = !r ^ + 1 . o^.) $ 

CD r r r V r / 


a 12 = 
cp 


22 


$ 

*z 

Ve 


r 

r 

r \r 

2 

r 2 


a£r]$cpq? 


"zCp 


<P 


rep 


x<- 

a af = 2$r$tcp “ 2aaqj 
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Again, as far the previously, on 2 it -will he more convenient 

to have the E(e) in terms of velocity components. Since again ^ = 0, 
they are 


E(z) 2w r w z ^ ®zz + {^ r w q?*r " (™r " &2 ) 3^3? $zcp + “ &2 ) "^T $zr + 

2 | w z^zz “ aa z^‘zz + x {^cp^zz^cpz + w z ^zcpj + v z$rz + w r^z$rz + 

~~2 j^cp “r^ $ qpcp“ aa z$cpcp] + r {^P^rz^rcp + w r + ' w r^rz^rr " ^z^rr j “ 


+ a *) 2 + a 2 $ rz - | w r 


] 


_(^p + oar) ^ + aaj (68a) 


<? 
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E(: 


■) = 2w r w z ^ $ zr + I Yr ' “ &2 ) 3?_|? $ cpr+ ^ " a2) S? *rr + 


2 |w z $ zr $ zz - aa r $ zz + i ^ zr <5 zC p+ ^ 


w 


$zcp “ 2w z — ^ 


+ cor 1 

S’zcpJ- 


*2 , . , i r ® c P r 

w r^zr + w z$£r^kr + 2 W( P 


3>c Dr ~ 2 ( w cp + «r) 


a 2 

®cpcp - aa^ ®cpcp+ — ®cpcp J + 


|[»<p4trV + w r V ~ 2 N * ^ «qj + T, r 4&. - aa^r} - 


1 V, ,.2,2 

+ mr) + a 


(68c) 

Finally, the derivatives of equations (68) to (70) with respect to v 
and of equation (68) with respect to u are required. 

In all these expressions, a and its first derivatives are obtained 

from 



*]- 


(r-i) i + 




V 


9‘ 


JL $ 2 $ 2 

r 2 v v r 


!) 
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CORCLUDIRG- REMARKS 

A method has been presented by means of which a detailed numerical 
description of many types of fluid flows in turbomachines can be obtained. 
The method is applied to one particular type described to the point 
where a solution is obtained merely by substitution of numerical values 
in the given equations . Once the algebra needed to obtain such equations 
has been done for any problem, the work involved in getting numerical 
results is not expected to be prohibitive, particularly in view of the 
fact that it largely consists of iterated integration, which type of 
calculation is well suited to rapid execution on high-speed computing 
machines . 


Lewis Flight Propulsion Laboratory 

Rational Advisory Committee for Aeronautics 
Cleveland, Ohio, February 12, 1952 
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APPENDIX - SYMBOLS 

The following symbols are used in this report: 

A,B,C,D constants in equation for m(z) 

fij' ^±i> 

Aj a • fourth-order minors obtained from the first four rows of 

,.r equation (10) 
a u a t 
*1 y *id 

A*3 coefficients in velocity transformation 

a sound velocity 


ry.2 a i3\ /*2 . a 2 


a 21 a 22 a 23 




'^-oxK 


v a 31 a 32 a 33 / \$ r $ z 


*£-« 

(^-Krf-e? 


r 2 >- 




$ $ 
r ^z 


id2-<or)* r 


r \r 


.2 2 

$r - a^ 


B 1 ^ 

B(€) 


coefficients in velocity transformation 

a 11 ^ + 2a 12 $> + 2a- 1 i3l> + a 22 fl> + 2a 23 I> + a 33 $ + b 

€ zz € zcp £ zr £ cpcp £ cpr £ rr £ 


1 

r 


» 2 ♦ a 2 




'12 


'34 


constants i = 1, • . . , 9 
intersection of 2-^ and 2g 
trailing edge 

_ 1 


p |_2 


21 + (cur)' 


nrv^*7 



2400 


NACA TN 2705 


37 


D± a coefficients in equation (39) 

* 

E(£) = -[I 33 ^ -2a 13 ]!: Q ze - [I 3 3p r - 2a 23 j^$ c p € +a 33 JL J r£ +B(t) 

£ = z^Cp, r 

F(|,tj) Legendre transformation of cp(z,r) in equation (49) 
f ^ integrands in solution of equivalent problem 

GjL coefficients in equation (39) 


g(z,r) 

given turning on £3 

H 

total enthalpy 

h 

static enthalpy 

I 

h + | W 2 - | o> 2 r 2 


_AS 

K 

e °V 

% 

constants, i = 1 , 2 

L(|) 

arbitrary function in equation (51) 

Z(z) 

intercept on cp = it /2 of generator of 2 

M 

Mach number 

m(z) 

angle with cp = 0 of generator of £ 

O 

N(^rj) 

integrand in equation (51) 

P 

static pressure 

r 

coordinate in cylindrical system 

3 

entropy 

s 

parameter on a surface 


T ' magnitude of tangent vector 

t parameter on a curve 
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u,v 


W 

v 

w z>V W r 


transformed coordinates 
•Wjp + cor 

magnitude of relative velocity 
parameter on a surface 
components of relative velocity 


z 

a 

P 

Ii 

r 

5i(z,r,v<p,W) 

€ 

e 


coordinate in cylindrical system 
angle between velocity vector and normal to £ 
complement of a, 
submatrices of D^ a 

ratio of specific beats 
coefficients in equation (59) 
an index denoting z, cp, or r 

independent variables of Legendre transformation 
Mach angle 

a function of entropy change across shock 


V- 

I 

Pi 



a root of equation (47c) 

roots of quadratic equation, i = 1, 2 


z 

Si 

22 

S 3 


initial surface 
shock surface 

suction surface ahead of shock surface 


suction surface behind shock surface 
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02 = (rqp z ) 2 + 1 + (r^) 2 

CT 1 = 02 " ( rC P z ) 2 ct| = a 2 -l a| = a 2 - (np r ) 2 


_2 

c 2 2 


a l = 

IV? ‘ K> 

- 

-2 



a 2 = 

1-P 2 


li 

nr 

Jp* - <^) 2 



roots of quadratic equation 


$ 

velocity potential 


cp 

coordinate in cylindrical system 


t 

function in equation (56) 


05 

angular velocity of wheel 



- ja 22 - 2a 32 qp r + a 33 cp 2 


Q 24 

a 11 - 2a“* r + a 33 zj 


1^14 

a 12 - a 1 ^. - a 23 ^ + a 3 ^^ 

, 

Subscripts: 


d 

downstream 


ij 

omitted columns in expansion of equation (10) 

H 

velocity component normal to shock 


TZ 

component tangent to shock in z = 

constant surface 

TR 

component tangent to shock in r = 

constant surface 
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1 ahead of shock surface 

2 behind shock surface 
Superscripts : 

ij distinguish coefficients in velocity transformations (52a) 

and (52b) 

(v) the approximation 

1 derivative with respect to t 
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